Newton's Laws of Motion in Form of Riccati Equation 



Marek Nowakowski and Haret C. Rosu 
Instituto de Ftsica de la Universidad de Guanajuato, Apdo Postal E-143, Leon, Guanajuato, Mexico 



O 
O 

a 

oo 



We discuss two applications of Riccati equation to Newton's laws of motion. The first one is the motion of a particle under 
the infiuence of a power law central potential V{r) = kr'^ . For zero total energy we show that the equation of motion can 
be cast in the Riccati form. We briefly show here an analogy to barotropic Friedmann-Robertson-Lemaitre cosmology where 
the expansion of the universe can be also shown to obey a Riccati equation. A second application in classical mechanics, 
where again the Riccati equation appears naturally, are problems involving quadratic friction. We use methods reminiscent to 
nonrelativistic supersymmetry to generalize and solve such problems. 

PACS number(s): 45.20.-d, 11.30 Pb 
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I. INTRODUCTION 

It is known that Riccati equations, in general of the 
type 



= f{x)y^ + 9{x)y + h{x) , 



(1) 



find surprisingly many applications in physics and math- 
ematics. For example, supersymmetric quantum me- 
chanics [Q, variational calculus ||^, nonlinear physics 
1^ , renormalization group equations for running coupling 
constants in quantum field theories Q| and thermody- 
namics 1^ are just a few topics where Riccati equations 
play a key role. The main reason for their ubiquity is 
that the change of function 



y 



dx 



(logz) 



(2) 



turns it into linear second-order differential equations of 
the form 



d z 



dx ^/dx 



1 dg 

2 dx 



dx 



log/ 



(3) 



that stand as basic mathematical background for many 
areas of physics. 

Since the Riccati equation is a widely studied nonlinear 
equation, knowing that the physical system under con- 
sideration can be brought into Riccati form has certainly 
many advantages. 

It is therefore of interest to look for yet different phys- 
ical problems which are governed by this first order non- 
linear equation. This can be a starting point to new av- 
enues in investigating analytical solutions of yet unsolved 
problems. In this paper we concentrate mainly on top- 
ics from classical mechanics and show that certain types 
of Newton's laws of motion are equivalent to the Riccati 
equation. 



II. THE POWER LAW CENTRAL POTENTIALS 

After implementation of the angular momentum con- 
servation law, the equation for the energy conservation 
in the case of a central potential V{r) is given by the 
standard expression 



h = —mr 
2 



/2 



2mr^ 



V{r) 



Taking a derivative with respect to time of Eq. 
into a second fundamental equation of the form 



(4) 
results 



dV{r) 
dr 



= 



Specializing from now on to a power law potential 
V{r) = kr^ , 



(5) 



(6) 



where k is the coupling constant and the exponent e can 
be either positive or negative, we obtain from (0) 



, , rnfr I 
V{r) = + 2 ■ 

Inserting the last equation in (0) gives 
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(7) 



(8) 



Under the assumption of — 0, this expression leads 
to a Riccati form. Obviously with E = 0, we restrict 
ourselves to the case k < 0. To explicitly derive from 
the Riccati equation we pass (as it is costummary in 
central potential problems) to an angle as a free variable 
(i.e., we consider r{9{t)). With 



dr 
d0 



and introducing 



r 
r 



(9) 



(10) 



1 



it can be readily shown, after some algebraic manipula- 
tions, that (ph reduces to 



(11) 



This is the Riccati equation for the motion of a particle 
in a central power law potential assuming E ^ 0. It is 
worth noting that no information about the coupling con- 
stant k enters the Riccati equation ( 11 ). Essentially what 
we have shown is that any solution of (Q) will also sat- 
isfy (0). The inverse is not necessarily true and should 
be examined in detail. Indeed, the coupling constant k 
should be explicitly contained in the solution for r{9) (see 
below). 

A special case which deserves to be briefly mentioned 
is e = —2. With this exponent, the choice i? = is, in 
general, only possible if ^ + fc < 0. Then directly from 

we conclude that ^ is a constant which, of course, is 
compatible with the Riccati equation (0). However, this 
constant cannot be determined by means of (^ij). This 
feature is also inherent in the general case. 

To discuss the case e 7^ 2, we first solve the Riccati 
equation (pTl). The solution can be easily found to be 



uj{9) = tan 



sin[(e + 2)6 + f3] 
cos[(e + 2)6 + 13] + I 



(12) 



where /3 plays a role of the integration constant. Going 
back to the definition of lu in (^|) we arrive at a solution 
for r{d) 



r{6) 



R 



[l + cos((e + 2)6» + /5)]i/(^+2) 



(13) 



where i? is a constant. As in the case e — —2 this constant 
can only be determined by inserting ( pi] ) into (^). The 
result is 



R = 



m\k\ 



(14) 



The last two equations represent then the analytical so- 
lution of the posed problem. We obtained this solution 
by transforming the original problem into a Riccati equa- 
tion. It might be that the laws of motion in Riccati form 
are only a curiosity. Given, however, the fact that only a 
few analytical solutions of the central potential problem 
are known, it is certainly a useful curiosity. Further- 
more, it is not excluded that this novel way opens new 
more general methods to solve problems in mechanics. In 
this context, we would like to mention here a yet differ- 
ent connection of the central potential problem with the 
Ermakov nonlinear differential equation H. We refer to 
the following form of the latter equation B 



which can be solved by the integrals 

dx r <f) 



(15) 



where a and b are integration constants and 
^(z) = 2 f f{z)dz . 



(16) 



(17) 



Taking p = const = m and suitably rescaling the distance 
r with the mass m, equation ( p^ is essentially identical 
to (^). Indeed, in this case the integrals in ( [1^ ) give 



t-U = — 



dri 



^2mV{ri) 



ii 

'■1 



(18) 



which, with a proper identification of 6, is the same as 
directly integrating (||). The interplay between the Er- 
makov equation and the central potential problems can 
be a useful tool of studying both problems. We con- 
jecture that certain invariants of the Ermakov equation 
could be also applied to the central potential problems. 



III. COSMOLOGICAL ANALOGY 

We want to point out here a beautiful but formal cos- 
mological analogy to the results of the previous section. 
We recall that in deriving the Riccati equation (9) we 
relied on a power law potential (^), a new parameter 6 
(the angle given in (^), and the assumption E = Q. The 
analogy to cosmology is based on these observations. In 
Friedmann-Robertson- Walker spacetime the set of Ein- 
stein's equations with the cosmological constant A set to 
zero reduce to differential equations for the scale factor 
a(t), which is a function of the comoving time t. Together 
with the conservation of energy-momentum tensor they 
are given by 

3S(t) = -47rG(p + Zp{p))a{t) (19) 
a(t)a(<) + 2d2(i) + 2k = 47rG(p - p{p))a^{t) (20) 

pa\t) = ^^{a\p+p{p)) . (21) 

In the above G is the Newtonian coupling constant, p is 
the pressure, p is the density and k can take the values 
0, ±1. Choosing the equation of state to be barotropic 



v{p) = (7 - i)p , 



(22) 



fixes essentially p to obey a power law behaviour of the 
form 



2 



(23) 



can be then deduced that the following equivalent form 
of 7 can be obtained 



and the remaining equation for a{t) reduce to a single 
equation, viz 



a{t) 



a{t) 



= 



c — —c — 1 
2 



(24) 



a{t)J a^{t) 
Introducing the conformal time 77 by 

dr] _ 1 

dt a{r/) ' 

it can be seen that ( |2^ ) is equivalent to a Riccati equation 

in the function u = where the dot means derivation 
with respect to 77 



(25) 



KC — 



(26) 



This cosmological Riccati equation has been previously 
obtained by Faraoni |^ and also discussed by Rosu in 
the context of late cosmological acceleration. The formal 
analogy to the mechanical case is obvious: the condition 
E = corresponds to A = 0, the angle 9 is replaced 
by the conformal time 77, and whereas in the mechanical 
example we had a power law behaviour of the potential, 
the barotropic equation of state forces upon p to satisfy 



p (X a 



-37 



As (nil) does not contain the coupling con- 



stant k, the cosmological Riccati equation (^6|) loses the 
information about G. 



IV. QUADRATIC FRICTION 

Starting with a constant force g (free fall, constant 
electric field, etc) and adding a quadratic friction with a 
positive friction coefficient > 0, we have, per excellence, 
a Riccati equation for the Newton's law of motion 

V ^ g ~ av^ , (27) 

with g = g /m and a = u /m. The general solution (which 
for reasons to be seen later in the text we denote by Vp) 
involves a free parameter A and reads 



Vp{t\g A) 



Ae- 



+ Ae-^ 



T = ag . (28) 



In the following we borrow some techniques from super- 
symmetric quantum mechanics. However, we do not fol- 
low strictly the supersymmetric scheme as the purposes 
in the quantum case and the mechanical case are quite 
different. We define a new time-dependent force by 

l{t;g ,a,A, Ai) =Vp{t;g ,a,X) +Xivl{t;g ,a,X) . (29) 



with a new parameter Ai > 0. We emphasize that (29) 



IS 

a definition given through the solution (28). From (27) it 



7 = 5 - {a- Xi)vl 



(30) 



This resembles supersymmetric quantum mechanics and 
we might be tempted to compare Vp to Witten's super- 
potential. 

To the new force we again add a quadratic function 
with a friction coefRcient Ai such that the new equation 
of motion becomes 



7 - '^1^' 



(31) 



This has the advantage that per construction Vp is a par- 
ticular solution of (|3l|). Equipped with this fact, one 
can proceed to construct the general solution which is 
a standard procedure in the general theory of the Ric- 
cati equation. Before doing so, it is instructive to dwell 
upon the physical meaning of the new force 7. Imposing 
g — (a — Ai)^ > 0, it can be seen that 7 > 0. Moreover, 
as obvious from (|2^ ) and (30), 7 goes to a constant for 
large t and has a kink-like behaviour. We can then envis- 
age a situation where 7 is a 'switch-on' function for a force 
becoming constant at some time. As mentioned above, 
by construction the problem (|^) is solvable because Vp 
is a particular solution of ( ^T] ) . By invoking the standard 
Bernoulli ansatz for the general solution Vg, namely 



1 



(32) 



we arrive at the differential equation (special case of the 
Bernoulli equation) for V 



Writing Vp as 



where 



V = 2Ait;pV" + Ai 



1 W„ 



Wp = e 



XiWp ' 



-Ai J Vpdt 



one is led to the solution for V 

Ai / W^dt + C 



V = 



The general solution is then given by 



''a = ''p+ X^Jwfdt + C 



(33) 



(34) 



(35) 



(36) 



(37) 



The initial value problem, w(0) = vq, is solved by fixing 
C through 



3 



Vo 



1-A 
1 + A 



1 

C 



(38) 



Up to integrals, the equation of motion (|31|) is solved. 
Setting A — e^^*^, we can rewrite ( p8| ) in the more conve- 
nient form 



Vp — — tanh(rt + 6) . 
a 

Then Wp can be computed explicitly 



[cosh(rt + (5)]^i/" 



(39) 



(40) 



It suffices to assume Ai = na, n G N leading to inte- 
grals of the type J cosli~"(x)(ia;, which can be solved in 
a closed analytical form by recursion formulae. Of course, 
the procedure outlined here can be generalized by start- 
ing with more complicated forces instead of the constant 
one. 
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IV. CONCLUSION 



In this paper we have pointed out the usefulness of the 
Riccati equation in studying certain mechanical prob- 
lems. We derived a Riccati equation for a central po- 
tential problem of the power law type assuming E = 0. 
This led us to an analytical solution of the problem. In 
a second step, we generalized the system of a constant 
force plus a quadratic friction to a time-dependent force 
and friction. We argued that this time-dependent force 
serves as a 'switch-on' function. The problem turned out 
to be solvable by means of a construction similar to su- 
persymmetric quantum mechanics. As indicated in the 
text, both applications can be generalized. 
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